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Summary. — We present a systematic study of the transverse-monientum depen- 
dent parton distributions (TMDs) in the framework of quark models. In particular, 
we review the general formalism for modeling the quark TMDs using the represen- 
tation in terms of overlap of nucleon light-cone wave functions. Such a formalism 
can be applied to a large class of quark models. We will discuss the building blocks 
of these different models, and we will use as explicit example for the calculation of 
the TMDs a light-cone constituent quark model. Within this model, we also pro- 
pose a phenomenological study of different observables, trying to learn how model 
parameters related to particular assumptions on the nucleon structure can be tuned 
to describe available experimental data. 
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1. Introduction 

The investigation of how the composite structure of a hadron, consisting of nearly 
niassless constituents, results from the underlying quark-gluon dynamics is a challenging 
problem, as it is of non-perturbativc nature, and displays many facets. What are the 
longitudinal momentum distributions of partons in a fast moving unpolarizcd or polar- 
ized hadron? What amount of transverse momentum do these partons carry, and how 
large is the resulting amount of orbital angular momentum? What is the spatial distri- 
bution of quarks inside a hadron as seen by a vector probe (coupling to the charge of the 
system), by an axial vector probe (coupling to the axial charge), or even seen by a more 
complicated probe? 
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Two new types of distributions such as the transverse-momentum dependent parton dis- 
tributions (TMDs) and the generalized parton distributions (GPDs) allow to address and 
quantify such questions. At leading-twist, the TMDs are defined as the diagonal matrix 
elements of the parton density matrix in the longitudinal and transverse-momentum 
space [l]-[5]. Taking into account all the independent spin-polarization states of the nu- 
cleon and partons, one obtains 8 independent TMDs describing the correlations in the 
three-dimensional momentum space between the nucleon and quark spin, and between 
the quark orbital angular momentum and both the nucleon and quark spin. On the other 
hand, GPDs are probability amplitudes related to the off-diagonal matrix elements of 
the parton density matrix in the longitudinal momentum space [6] -[17]. After a Fourier 
transform to the impact-parameter space, they also provide a three-dimensional picture 
of the hadron in a mixed momentum-coordinate space [18]-[20]. 

The information encoded in the TMDs and the GPDs can be collected in a unified 
framework through the generalized transverse-momentum dependent parton distribu- 
tions (GTMDs); for a classification see Refs. [21, 22]. The GTMDs contain the most 
general one-body information about partons, corresponding to the full one-quark density 
matrix in momentum space [23] . GTMDs depend on the three-momentum of the quarks 
as well as on the four-momentum A which is transferred by the probe to the nucleon. 
The Fourier transform of the GTMDs to the impact-parameter space allows to define 
the Wigner distributions of the parton-hadron system, which represent the quantum- 
mechanical analogues of the classical phase-space distributions [24, 25, 26, 16]. Wigner 
distributions provide 5-dimensional (two position and three momentum coordinates) im- 
ages of the nucleon as seen in the infinite-momentum frame. As such they contain the full 
correlations between the quark transverse position and three-momentum. In particular 
limits, they reduce to different parton distributions and form factors as shown in Fig. 1. 
The different arrows in this figure represent particular projections in hadron and quark 
momentum space, and give the links between the matrix elements of different reduced 
matrices. In the forward limit A = 0, they reduce to the TMDs, while the integration 
over the quark transverse momentum fcj^ leads to the GPDs. The common limit of TMDs 
and GPDs is given by the standard parton distribution functions (PDFs). The integra- 
tion over the longitudinal-momentum fraction x brings to the lower plane of the box in 
Fig. 1, and at each vertex we find the restricted version of the operator defining the dis- 
tributions in the upper plane. For example, in correspondence of the TMDs we encounter 
the transverse-momentum dependent spin densities (TMSD), while from the GPDs we 
obtain the form factors (FFs). Both FFs and TMSDs have the charges as common limit. 
Each of these observable is sensitive to specific properties of the nucleon structure, and 
contributes from different perspective to form a global and comprehensive view [27]. In 
these lectures we will concentrate on the model calculation of TMDs, restricting our- 
selves to the quark contribution. We will exploit the relations with different observables 
and parton distributions in the box of Fig. 1 to test different model assumptions and 
their implications for the partonic structure of the nucleon. A variety low-energy QCD- 
inspired models have been employed to model the TMDs [23], [28]-[55]. Although they 
all have in common that they strongly oversimplify the complexity of the QCD dynamics 
in hadrons, studies in different models based on often complementary assumptions, help 
to unravel non-perturbative aspects of TMDs. Insights into non-perturbative proper- 
ties are of particular interest when confirmed in various models. The practical value of 
model results is that they can be used to predict new observables, or to guide educated 
Ansatze for fits of TMD parametrizations. Especially in the context of TMDs one should 
not underestimate the conceptual importance of model calculations. Model calculations 
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Fig. 1. - Representation of the projections of the GTMDs into parton distributions and form 
factors. The arrows correspond to different reductions in the hadron and quark momentum 
space: the solid (red) arrows give the forward limit in the hadron momentum, the dotted 
(black) arrows correspond to integrating over the quark transverse momentum and the dashed 
(blue) arrows project out the longitudinal momentum of quarks. 



demonstrated the existence of effects [56], paved the way towards an understanding of 
universality in the fragmentation process [57], established new TMDs [58, 59], see [60] 
for a review. 

The manuscript is organized as follows. In sect. 2 we review the definitions of the 
leading-twist TMDs and introduce a convenient representation of the quark-quark cor- 
relator in terms of the net-polarization states of the quark and the hadron. In sect. 3 
we discuss model relations among TMDs which hold in a large class of quark models. In 
particular, we review the derivation given in Ref. [61] to explain the physical origin of 
such relations. There are two linear relations and a quadratic relation which are flavour 
independent and involve polarized TMDs, while a further linear relation is flavour depen- 
dent and involves both polarized and unpolarized TMDs. The relations among polarized 
TMDs connect the distributions of quarks inside the nucleon for different configurations 
of the polarization states of the hadron and the parton. As a consequence, it is natural 
to expect that they can originate from rotational invariance of the polarization states of 
the system. Rotations are more easily discussed in the basis of canonical spin. Therefore, 
instead of working in the standard basis of light-cone helicity, we introduce in sects. 4'1 
and 4'2 the tensor correlator defining the TMDs in the canonical-spin basis. Such a 
representation is used in sect. 4'3 to discuss the consequences of rotational symmetries of 
the system. In such a way we will be able to identify the key ingredients for the existence 
of relations among polarized TMDs in quark models. 

In order to complete the discussion, including the flavour-dependent relation among 
polarized and unpolarized TMDs, we need to introduce specific assumptions about the 
spin-isospin structure of the nucleon state. Therefore, in sect. 5 we discuss the conse- 
quences of rotational invariance using the explicit representation of the TMDs in terms 
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of three-quark (3Q) wave functions. In particular, in sect. 5'1 we derive the overlap 
representation of the TMDs in terms of light-cone wave functions (LCWFs), while the 
corresponding representation in terms of wave function in the canonical-spin basis is 
given in sect. 5'2. In sect. 5'3 we discuss the constraints of rotational symmetry on the 
nucleon wave function and, as a result, we give an alternative derivation of the flavour- 
independent relations among TMDs. Finally, in sect. 5'4 we discuss the constraint due to 
SU{6) symmetry of the spin-flavour dependent part of the nucleon wave function. This 
additional ingredient allows us to explain the origin of the flavour-dependent relation. 
The formal derivation of the relations among TMDs is made explicit within different 
quark models in sect. 6. There we review different quark models which have been used 
in the literature for the calculation of TMDs. In particular, we discuss the key ingredients 
of the models, showing how and to which extent the conditions which lead to relations 
among TMDs are realized. In sect. 7 we also study the phenomenological implications 
of the breaking of SU{6) symmetry. We will work in the framework of a light-cone con- 
stituent quark model, giving the essential step for the explicit calculation of the TMDs 
in the model in sect. 7'1. In sect. 7'2 we discuss the model results for T-even TMDs 
and related spin densities in the transverse-momentum space. In sect. 7'3 wc propose 
a phenomenological study of different observables, fitting the model parameter related 
to iS'C/(6)-symmetry breaking to the double spin asymmetry in inclusive deep inelastic 
scattering (DIS). The results from this fit are then exploited to obtain predictions for 
other observables which are particularly sensitive to S'C/(6)-breaking effects, like the dou- 
ble spin asymmetry for semi-inclusive deep inelastic scattering (SIDIS) and the neutron 
electric form factor. A summary of our findings is given in the final section. Techni- 
cal details and further explanations about the representation in terms of nucleon wave 
functions are collected in three appendices. 

2. Transverse-momentum dependent parton distributions 

2'1. Definitions. - In this section, we review the formalism for the definition of TMDs, 
following the conventions of refs. [4, 62, 63, 64]. Introducing two lightlike four- vectors 
n± satisfying n+ • n_ = 1, we write the light-cone components of a generic four- vector a 
as [a'^ , a~ , a±] with a* = a-Uzp. The density of quarks can be defined from the following 
quark-quark correlator (^) 
(1) 



(27r) 



where /c+ = xP^, tjj is the quark field operator with a, 6 indices in the Dirac space, 
and U is the Wilson line which ensures colour gauge invariance [65]. The target state 
is characterized by its four-momentum P and covariant spin four-vector S satisfying 
P^ = M^, S^ — —1, and P ■ S ^ 0. We choose a frame where the hadron momentum has 



no transverse components, P 



and so S ~ 



c P^ _q M a 

^ M ' ^ 2P+ ' -' 



with 



S = 1 . From now on, we replace the dependence on the covariant spin four- vector S by 
the dependence on the unit three- vector S = (S'_l, Sz)- 



(^) In order to simplify the notation, we will omit the colour index, keeping in mind that the 
quark-quark correlator is diagonal in the colour space. 
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TMDs enter the general Lorentz-covariant decomposition of the correlator ^ab{x, k±, S) 
which, at twist-two level and for a spin-1/2 target, reads 

(2) ^x, k^, S) = i{/i fi+ - ^i^i^ ftr ii+ + S, 9iL 75?i+ + ^^ 9it 75?^+ 

where e^^ = ^^t^ = 1j ^-i^d the transverse four-vectors arc defined as a^_ = [0,0, aj^]. 
The nomenclature of the distribution functions follows closely that of Ref . [4] , sometimes 
referred to as "Amsterdam notation" : / refers to unpolarized quarks; g and h to lon- 
gitudinally and transversely polarized quark, respectively; a subscript 1 is given to the 
twist-two functions; subscripts L or T refer to the connection with the hadron spin being 
longitudinal or transverse; and a symbol _L signals the explicit presence of transverse 
momenta with an uncontracted index. Among these eight distributions, the so-called 
Boer-Mulders function h^ [62] and Sivers function f^ [66] are T-odd, i.e. they change 
sign under "naive time-reversal" , which is defined as usual time-reversal but without in- 
terchange of initial and final states. All the TMDs depend on x and k\. These functions 
can be individually isolated by performing traces of the correlator with suitable Dirac 
matrices. Using the abbreviation ^I'"' s Tr[$r]/2, wc have 

(3) $[^+l(a;, fci, S) = h- '-^j^ /iV, 

(4) ^lf^i^][:,^k^,S) = S,g^L + ^^9iT. 

(5) $[^-'^''^1 (x, k^,S) = Si h, + S,^ hi^ + 51 ^'^^%-f^'" hij, + '^hi, 

where j = 1, 2 is a transverse index, and 

(6) hi = hiT + 2]v7T h^rj.. 

The correlation function $'''' '(x, fc_L,S^) is just the unpolarized quark distribution, 
which integrated over fc_L gives the familiar light-cone momentum distribution fi{x). All 
the other TMDs characterize the strength of different spin-spin and spin-orbit correla- 
tions. The precise form of this correlation is given by the prefactors of the TMDs in 
Eqs. (3)- (5). In particular, the TMDs qil and hi describe the strength of a correlation 
between a longitudinal/transverse target polarization and a longitudinal/transverse par- 
ton polarization. After integration over k±^ , they reduce to the helicity and transversity 
distributions, respectively. By definition, the spin-orbit correlations described by f^, 
giT, hi , h^^ and h^rp involve the transverse parton momentum and the polarization of 
both the parton and the target, and vanish upon integration over fcj_. 

If one calculates these distributions in the light-cone gauge A+ = using advanced 
boundary condition for the transverse component of the gauge field, the gauge links in the 
quark-quark correlator can be ignored [67]- [69]. However, in this case the wave function 
amplitudes are not real and, apart from the structural information on the hadron, they 
also have an imaginary phase mimicking the final state interactions [70, 71]. 
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2'2. Helicity and four- component bases. - The physical meaning of the correlations 
encoded in TMDs becomes especially transparent when using for the quark field operator 
the Fourier expansion in momentum space in terms of light-cone Fock operators. To this 
aim, we first introduce the Fourier expansion in momentum space of the quark field 
operator ■0' of flavour q (see emphe.g. [72]) 

— — -j7;^e(fc+) 2_]{qA(fc)wLc(^, A) exp(-ifc+z" + ik_L ■ z_l) 
"" A 

(7) +ql{k)vLc{k) exp(+ifc+z" - i'kj_ ■ z_l)}, 

where k — (fc''',kj_), A, represents the parton helicity, and q and q' are the annihilation 
operator of the quark fields and the creation operators of the antiquark fields, respectively, 
which fulfill the anticommutations relations 

(8) {qyCk'),qlCk)} ^ {qx'Ck'),qiCk)} = lQ^'k^S{k'+ - k+)S^^\k\ -k^)Sg,,Syx. 

In Eq. (7), the quark LC spinors ULcik,X) and VLc{k,\) are the free Dirac light-cone 
spinor of the quark and antiquark, respectively (see Appendix A) . 

For simplicity, in this study we will consider only the quark contribution to TMDs, 
ignoring the contribution from gauge fields and therefore reducing the gauge links in 
Eq. (1) to the identity. As a consequence, we will not be able to discuss the T-odd TMDs, 
which vanish identically in absence of gauge-field degrees of freedom. Furthermore, we 
decompose the correlator (1) into the different quark fiavour contributions 

(9) * = E*«- 

Following the lines of refs. [73]- [75], we obtain at twist-two level 

(10) $fl(x,fc^,5) = ^(P,Sl^gl,(fc)qA(fc)MFl^'^|P,S'), 

A' A 

where N ^ [2a;(27r)3]^ j(3)(0) and I/FJa'a ^ u^c:(fc, A')ruLc(fc, A)/2fc+. As a resuh, 
using (3)-(5), the quark contribution to the T-even TMDs is given by 

(11) n{x,kl)^{P,K\V\k)\P,K), 

(12) A5?^(.T,fci) = (P,AlA«(fc)lP,A), 

(13) ^i_^g^^(:c,fci) = {P,SM'mP.S^), 

(14) A^;i^^«(x,fci) = i^sign(A)(F,Alri(fc)|P,A), 

A 

(15) S^/^?^(a;,fci) + ^i^^/.^^'(x,fcl) = i^sign(5^)(P,5j.lTi(fc)|P,5^). 
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where the quark operators V^, A"^, and T"^ arc defined as 

(16) y^(fc) = ^9l(%A(fc), 

A 

(17) A^k) = Y.'gniX)qlCk)q,ik), 

A 

(18) T^^Ck)=[Tl{k)]^ ^2qlCk)q-Ck)- 



The quark operators V and A"^ have a probabihstic interpretation since they are writ- 
ten just in terms of number operators 7V^ = ql^(k)qx{k). The operator T^ has also 
a probabilistic interpretation but only when written in terms of transversely polarized 
operators 



(19) 



s± ■ Tl = ^sgn(si)glj^)g, Jfc), 



where for a generic direction s = (sin 0s cos (f>s, sin^^ sinc/)^, cos 0s), the creation operators 
for a polarized quark can be expressed in terms of creation operators for a longitudinally 
polarized quark as 

(20) {qU qU) = {ql ql-)^iOs,<f>s), 

where the SU{2) rotation matrix u{9, (p) is given by 



(21) 



(a A\ i cos f e-*'^/^ 
"(^''^)=( sinfe^^/2 






The operators V^k), A''{k), and s±_ ■ T\^{k) give the number, the effective longitudinal 
polarization, and the effective transverse polarization of quarks with flavour q and light- 
cone momentum fc, respectively. 

We can further disentangle (15) by choosing appropriate target and quark polariza- 
tions 



(22) 
hl{x,k\) 



J2 sign iS.KP, S,\T^Mp, 5,) + ^ sign {Sy){P, Sy\T^{k)\P, Sy) 



(23) 



M2 



— hj^rp{x,kj_) — - 



J2 sign {S.){P, S.\n{k)\P, S,) - Y. sign {Sy){P, Sy\T^{k)\P, Sy 



or 



(24) 



hi:S{x,kl) ^lY.^ign{S,){P,S,\T^{k)\P,S,) = lj2^ign{Sy){P,Sy\TSCk)\P,Sy), 
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where the nucleon polarization state \P, S) in a generic direction S = (sin ^5 cos 05, sin Os sin (ps, cos Os) 
can be written in terms of longitudinally polarized states |-P, A) as 



(25) 



(|P,+5), |P,-5)) = (|P,+>, |P,-))«(%,, 



In the literature, one often represents correlators in terms of helicity amplitudes which 
treat in a symmetric way both quark and target polarizations 



(26) 



$ 



1 



'\J 



l,y^^^{x,k^) = j^{P,A'\ql,ik)q^ik)\P,A) 



Thanks to light-cone parity invariance, the helicity amplitudes are not all independent 
and are constrained by the following relation 



(27) 



^I'X'.AxCk) - '^'-A'-X'-A-xikp), 



where kp = (x, — fc^, ky). From the definitions in Eqs. (11)-(15) it follows that it is pos- 
sible to express each TMD as linear combination of helicity amplitudes as described in 
Table I. In particular we note that each TMD involves combinations of helicity ampli- 
tudes with a definite value of orbital angular momentum transfer A?^ between the initial 
and final state. This is a consequence of the conservation of the total angular momen- 
tum which gives the constraint A/^ = (AA — AA)/2, with AA = A' — A and analogously 
AA = A' - A. 

Table I. - Relations between the 8 leading-twist TMDs and helicity amplitudes. The table can be 
used to obtain the TMDs as functions of the helicity amplitudes as well as the inverse relation. 
The entries in the table give the coefficients in the linear combination which relates the two sets 
in the first column and in the first row. The helicity amplitudes are also classified in terms of 
the transfer of orbital angular momentum A/2 between the initial and final state. 



n+,++ 






$! 



$! 



A/2 = +1 

$9^ 



A/2 = +2 



n 

hi 

r±q 
J IT 


1 
1 


q 

9iT 




ht' 




htl 




"it 





1 

-1 



■ Mkf 



k 

Mk 
k 



R 



J_ 



■ Mkf, 



Mkn 



~kr 



{"^y 



We can collect the information encoded in Table I by introducing the following matrix 
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form 

(28) 



^A'A'.AA 










( Hf!+9h) 


-^{^h^'-hil) 


jt{^flT' + 9lT) 


hi 




tr {^hi^ + hil) 


kJf?-9L) 


2KP "it 


tl (*/it' - 91t) 




2M (*/lT 9it) 


k1 u-i-q 

2A/2 IT-IT 


hin-9iL) 


-^{^hi^ + hil 




[ K 


-|^(^/iV^ + 3?t) 


^ (*'h " - hii) 


hiri+9lL) 



where the entries in the rows correspond to (A'A') = (++), (H — ), ( — h), ( ) and the 

entries in the eolumn are hkewise (AA) = (++), (H — ), ( — h), ( )• 

Alternatively, we can represent the quark correlator (10) in the four-component basis. 
There are only four twist-two Dirac structures rtwist-2 = {7^7 icr^^75: icr^^Ts: 7^75}, 
corresponding to the four kinds of transition the light-cone helicity of the active quark 
can undergo (see e.g. [76]-[78]) 

(29) Afl-'^l^'^ = d^'^, A/[»'^^+75]A'A ^ (ajf^, A/[t+75]a'a ^ ^^.^^a'a 

with CTi the three Pauli matrices. For further convenience, we associate a four-component 
vector (^) to every quantity with superscript F 

(30) ali^l ^ a'' = {a'^.a^a^a') = (al^^l, a'^'^'^^^U'*"'^''^!, a'^^^^l) . 
With this notation, the correspondence (29) takes the simple form 

(31) M"^'^ ^ i^'")^'^, 

where a'^ = (lj(?)- The symbols ct^ and cr^ satisfy the relations 

(32) 1 ian''\cT,)rr' = S^:S^, ixr [aV,] = i ^(a'')'''K)AA' = -5^ 

A'A 

One can think of (j''^ ^ = (a-'')^ ^ as the matrix of a mere change of basis, the one being 
labeled by the couple AA' and the other by i' 

(33) a^'-^a'^^'^aAA'. 

A'A 

In this basis, we can introduce the tensor correlator ^^'^ {x, k±) which is related to helicity 
amplitudes as follows 

(34) *r = ^ E (-")^^' (-'')''' *1'A',AA: *I'A',AA = ^<J'rK)A'AK)AA'. 

A' A A' A 

(^) Note this is not a Lorentz four- vector but Einstein's summation convention still applies. 
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The tensor correlator is then given by the following combinations of TMDs 



g 



(35) 



ky_ ri-q 
M -I IT 

_fcx f-L« 
M JIT 

\ 

( n 

fcy f^q 
M J IT 



' M JlT 





hi 



ky_ ,±q 
kl-kl 



2M^ "it 



"-IT 



M "-IL 

M 'n 






/.? 



k: 



"'IT 



"IT ""a; 

k^uM 
M 'hL 



""IT y 



it) ^ y 



2M^ "it 
A/ "iL 



^ 
fcx „9 

fcy 9 

M 9lT 

9\l > 



M "l 

fcy , _Lg 
M ''■IL 



\ 



fcx 9 

M yiT 

fcy <? 
M ^IT 

5?^ / 



where we introduced the notations h^^ = h\± j^ h^^ and fc^ = ki/kj_ with fcj^ = |fc±l- 
The component $°° gives the density of quarks in the target irrespective of any polariza- 
tion, i.e. the density of unpolarized quarks in the unpolarized target. The components 
<f>^^ give the net density of quarks with light-cone polarization in the direction Ej in the 
unpolarized target, while the components <i>5' give the net density of unpolarized quarks 
in the target with light-cone polarization in the direction e^. Finally, the components 
$*■' give the net density of quarks with light-cone polarization in the direction Bj in the 
target with light-cone polarization in the direction e^. The density of quarks with def- 
inite light-cone polarization in the direction s inside the target with definite light-cone 
polarization in the direction S is then obviously given by $g(x, fej^, S, s) = -^ S^^^'^s^, 
where we have introduced the four-component vectors 5^ = (1, S) and Si, = (1, s). 

3. Model relations 

In QCD, the eight TMDs are all independent. It appeared however in a large panel of 
low-energy models that relations among some TMDs exist. At twist- two level, there are 
three flavour-independent relations (^), two are linear and one is quadratic in the TMDs 



(36) 

(37) 
(38) 



5?L 



fflT 



(5?t) 



"■1 ^ 2A/2 "it 

hi! = 0, 

^ ' 2hl hl^ = 



= 0, 



A further flavour-dependent relation involves both polarized and unpolarized TMDs 

(39) V'ifl + glp = 2hl 

where, for a proton target, the flavour factors with q = u,d are given by P" = | and 
T^'^ = —k- As discussed in Ref. [48], at variance with the relations (36)- (38), the flavour 
dependence in the relation (39) requires specific assumptions for the spin-isospin structure 
of the nucleon state, like SU{6) spin-flavour symmetry. 



{^) Other expressions can be found in the literature, but are just combinations of tlie relations 
(36)-(38). 
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A discussion on how general these relations are can be found in Ref. [48]. Let us 
just mention that they were observed in the bag model [48]-[41], light-cone constituent 
quark models [75], some quark-diquark models [28]- [39], the covariant parton model [37] 
and more recently in the light-cone version of the chiral quark-soliton model [23] . Note 
however that there also exist models where the relations are not satisfied, like in some 
versions of the spectator model [47] and the quark-target model [51]. 

As already emphasized, the model relations (36)-(39) are not expected to hold iden- 
tically in QCD, since the TMDs in these relations follow different evolution patterns. 
This implies that even if the relations are satisfied at some (low) scale, they would not 
hold anymore for other (higher) scales. The interest in these relations is therefore purely 
phenomenological. Experiments provide more and more data on observables related to 
TMDs, but need inputs from educated models and parametrizations for the extraction 
of these distributions. It is therefore particularly interesting to see to what extent the 
relations (36)-(39) can be useful as approximate relations, which provide simplified and 
intuitive notions for the interpretation of the data. Note that some preliminary calcula- 
tions in lattice QCD give indications that the relation (37) may indeed be approximately 
satisfied [79]-[81]. 

Using two different approaches, we show in the next sections that the flavour-independent 
relations (36)-(38) can easily be derived, once the following assumptions are made: 

1. the probed quark behaves as if it does not interact directly with the other partons 
(i.e. one works within the standard impulse approximation) and there are no 
explicit gluons; 

2. the quark light-cone and canonical polarizations are related by a rotation with axis 
orthogonal to both light-cone and k± directions; 

3. the target has spherical symmetry in the canonical-spin basis. 

From these assumptions, one realizes that the flavour- independent relations have essen- 
tially a geometrical origin, as was already suspected in the context of the bag model 
almost a decade ago [82]. We note however that the spherical symmetry is a sufficient 
but not necessary condition for the validity of the individual flavour-independent rela- 
tions. As discussed in the following section, a subset of relations can be derived using 
less restrictive conditions, like axial symmetry about a specific direction. 

For the fiavour-depcndent relation (39), we need a further condition for the spin- 
flavour dependent part of the nucleon wave function. Spcciflcally, we require 

4. SU{6) spin-flavour symmetry of the wave function. 

As shown in sect. 6, it is found that all the models satisfying the relations also satisfy 
the above conditions. We are not aware of any model calculation which satisfies some or 
all the three flavour-independent relations and at the same time breaks at least one of 
the conditions 1-3. However, this is not a priori excluded. 

4. Amplitude approach 

The first derivation of the TMD relations stays at the level of the amplitudes. As we 
have seen, the TMDs can be expressed in simple terms using light-cone polarization. On 
the other hand, rotational symmetry is easier to handle in terms of canonical polarization, 
which is the natural one in the instant form. We therefore write the TMDs in the 
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canonical-spin basis, and then impose spherical symmetry. But before that, we need to 
know how to connect light-cone helicity to canonical spin. 

4'1. Connection between light-cone helicity and canonical spin. - Relating in general 
light-cone helicity with canonical spin is usually quite complicated, as the dynamics is 
involved. Fortunately, the common approach in quark models is to assume that the target 
can be described by quarks without mutual interactions. In this case the connection 
simply reduces to a rotation in polarization space with axis orthogonal to both k± and 
Ez directions. The quark creation operator with canonical spin a can then be written in 
terms of quark creation operators with light-cone helicity A as follows 



(40) 



= E^. 



(1/2)* 
A 



^1 



with 



D 



(1/2)* _ 



fcL sin I 



kr sm I 
cos I 



Note that the rotation does not depend on the quark flavour. The angle 6 between 
light-cone and canonical polarizations is usually a complicated function of the quark 
momentum k and is specific to each model. It contains part of the model dynamics. The 
only general property is that — > as fcj^ — > 0. Due to our choice of reference frame 
where the target has no transverse momentum, the light-cone helicity and canonical spin 
of the target can be identified, at variance with the quark polarizations. 

4'2. TMDs in canonical- spin basis. - The four-component notation introduced in 
sect. 2'2 is very convenient for discussing the rotation between canonical spin and light- 
cone helicity at the amplitude level. Since the light-cone helicity and canonical spin of 
the target can be identified in our choice of reference frame, we expect the canonical 
tensor correlator ^q' to be related to the light-cone one in Eq. (35) as follows 



(41) 



Cq q ^p 



with O some orthogonal matrix O^O ~ 1 representing the rotation at the amplitude 
level. From Eqs. (26), (34), (32) and (40) we find that the orthogonal matrix is given by 



P 2 L '^ 

^ ^^ r)(i/2) / \ 7-,vv^;* t^yy 



,(1/2)* 



2 ^^ 

(t'ctA'A 



(42) 



/I 






kl + kl cos 9 




—kxky (1 — COS 6*) 


\ 

~kx sin 9 





— kxky (1 — cos9) 


kl + kl cos 9 


^ky sin 9 


\o 


kx sa\9 


ky sin 9 


cos 9 ) 



The canonical tensor correlator then takes the form 



(43) 





( f! 


ky 1,-i-q 


A/ 1 


\ 


^Cq - 


ky rJ-q 
M JIT 


Hit ""x "T '^iT y 


(vlT ~ IT) kxky 


k^<i 
M Bit 


-tf^T' 


(Oir ~ '*it) ^xky 


"it "^x ^ ^Ht S 


M Bit 




\ 


M ^IL 


kjLU-^q 
M ^IL 


0?L / 
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where we introduced the notations 



(44) 



(45) 



M ^IL 


M 01T 
UlT 



/cos 9 
l^sin6' 

/cos 9 

l^sin e 



cost 



cost 



91l 

M "-IL 

M yiT 

"-IT 



Comparing Eq. (43) with Eq. (35), wc observe that the multipole structure is conserved 
under the rotation (41). The rotation from hght-cone to canonical polarizations affects 
only some of the multipole magnitudes, see Eqs. (44) and (45). 

Note that the orientation of the axes in the transverse plane has been fixed arbitrarily. 
There is however a privileged direction given by the active quark transverse momentum 
fe j^ . Choosing the orientation of transverse axes so that cither k±^ ~ kj_ e^ or k± = 
k± By simplifies the transformation, as it eliminates the cumbersome factors kx and ky 
in Eqs. (42) and (43). Choosing e.g the second option, the orthogonal matrix of Ec^. (42) 
reduces to 



(46) 



On 



/I 




Vo 





1 
cos( 
sin ^ 



\ 


— sin I 

cos 6' / 



and the light-cone and canonical tensor correlators take the following simpler forms 



(47) 



(48) 



9 



'^'J, 



( n 

k± fJ-q 
M -I IT 



V 
/ f? 

k± f-i-q 
M -I IT 



V 



M "l 
"■IT 





ki_ i,-L9 
A/ "l 

"■IT 






"■IT 






^ 


M yiT 
9\l ' 









+q 

IT 


k± „q 

0?L 



Playing a little bit with Eqs. (44) and (45), we find 



(49) 




/cos 9 — sin 
lsin0 cos 




M \3lT + ^IL ^ 
01L MiT 



and three expressions invariant under the rotation (41) 

(50) (^5?T)' + (/^S)'-(^fl?T)' + (b^ 

(51) i9lLf+{^htiy = i9l,f+{!^hi 



(52) 



9il9it + "iL "iT — 01L 01T + Oil ■'it- 



These three invariant expressions have a simple geometric interpretation. The three- 
component vector ^^ S''^*^ = tt^ represents the net light-cone polarization of a quark 
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with three- momentuni {xP~^ , k±) and flavour g in a target with net polarization in the 
S'-dircction. From Eq. (41), we see that the vector Tr^g representing the net canonical 
polarization of the quark is obtained by a rotation of tt^ 



(53) 



'cs 



E'^: 



|fcQfcj; 



It follows automatically that tt^ • tt^, is invariant under the rotation (41) 



(54) 



r« - 



'cs 



'cs'- 



Expressions (50) and (51) are obtained from (54) for the cases S ^ S' = e± and S = 
S = e^, respectively. They just express the fact that the magnitude of quark net 
polarization is the same in the light-cone helicity and canonical-spin bases. Expression 
(52) is obtained for the case 5' = ej^ and 5'' = Bz- All the remaining cases do not lead 
to new independent expressions. 



4'3. Spherical symmetry. - We are now ready to discuss the implications of spherical 
symmetry in the canonical-spin basis. Spherical symmetry means that the canonical 
tensor correlator has to be invariant OJ^^cqOr — $cg under any spatial rotation On = 
(or) ■^ith R the ordinary 3x3 rotation matrix. It is equivalent to the statement 
that the tensor correlator has to commute with all the elements of the rotation group 
^CqOR ~ Oji^cq- As a result of Schur's lemma, the canonical tensor correlator must 
have the following structure 



(55) 



*&; 

























By 



Comparing this with Eqs. (43) or (48), we conclude that spherical symmetry implies 



(56) 
(57) 
(58) 



n = A^, 



^' -t)t^^h 



ilL 



'IT 



a'' 



Kl = f 



IT 



9 - D? 
IT ~ ^ ' 



0. 



Clearly, only the monopole structures in the canonical-spin basis arc allowed to survive. 
Furthermore, the Sivers and Boer-Mulders functions f^ji and h^ '^ vanish identically, as 
expected from the fact that we are neglecting gauge-field degrees of freedom. 

Note however that the monopole structures in the canonical-spin basis generate higher 
multipole structures in the light-cone helicity basis. It follows that spherical symmetry 
imposes some relations among the multipole structures in the light-cone helicity basis, 
and therefore among the TMDs. Inserting the constraints (57) and (58) into Eq. (49), 
we automatically obtain the linear relations 



(59) 



'k^ l^t ^ i^U 



M \9iT 



IL _ 



"■IT 
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Using now the constraints from spherical symmetry in Eq. (50), we obtain the quadratic 
relation (38) 

(60) = (^5?t)' + {K^f {h-^f - ft [(.9?t)' + ^hlht^ 

The linear relations (36) and (37) being satisfied, the Eqs. (51) and (52) do not lead to 
independent quadratic relations. 

We have seen that spherical symmetry is a sufficient condition (^) to obtain all three 
flavour-independent relations. Restricting ourselves to axial symmetries, we find that 
some of the relations can already be obtained. For example, axial symmetry about e^ 
alone implies the quadratic relation (38) and 

(61) gl^gl^ + hilhtS = 0. 

Axial symmetry about k± x e^ implies the two linear relations (36) and (37). The relation 

(61) is naturally also satisfied but is not independent. 

5. Wave- function approach 

Many quark models are based on a wave- function approach. We therefore translate 
here the derivation of the previous section in the language of 3Q wave functions. The 
advantage is that we can then also discuss the additional SU{6) spin-flavour symmetry 
needed for the flavour-dependent relation (39). 

5'1. Overlap representation of the TMDs on the light cone. ~ Restricting ourselves to 
the 3Q Fock sector, the target state with definite four- momentum P = [P+,^p^,Oj^] 
and light-cone helicity A can be written as follows 

(62) |P,A)= J2 E /[dx]3[d2fc^]3^A;A/r(fci,fc2,^3)|{A.,g.,feO), 

A1A2A3 919293 

where V'a'a^'a'^(^1i^2, ^3) is the three-quark light-cone wave function (3Q LCWF) with 
Xi , Qi and ki referring to the light-cone helicity, flavour and light-cone momentum of quark 
i, respectively. The total orbital angular momentum of a given component ip^ _,j^ j,,^ is 
given by the expression l^ = A — Ai — A2 — A3 with A, Aj = ±i. The integration measures 



(*) From Eqs. (49) and (50), one can see that the minimal conditions are actually 

0?L - 1)S = 0, 



,9 j_ h^l 

^2 ,,-„.2 



ilT ^ ^IL — 0' 



['rrslT) +{^tS) -{h;S) =0. 

They are indeed fulfilled by spherical symmetry, see Eqs. (57) and (58). 
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in Eq. (62) are defined as 

[dxU = 



(63) 



" 3 

[[dx. 


/ 
S 
\ 


1 r d ki± 



2(27r)=^<5(2)K]fc,^ 



^i=l 



Choosing to label the active quark with i = 1, the TMDs can be obtained by the following 
overlaps {^) of 3Q LCWFs 



(64a) 
(64b) 
(64c) 
(64d) 
(64e) 
(64f) 
(64g) 
(64h) 



f! 



91l 



hi 



M J IT 



M 'JlT 



d[23 
d[23 
d[23 
d[23 
d[23 



'^hi'^= I d[23 



^hil= /d[23: 



SIP" "it 



d[23 



EE[\<ZT\" + \^XT\"] 



•^2-^3 92 93 



EEo^:a"i'-i^-Axn 



A2A3 9293 



A2A3 9293 



99293 
A2A3 ' 



A2A3 <3'2Q'3 



QQ2q3 
+ A2A3 



EE25ie[fci.(^:i,,3 
A2A3 9293 



992 93 V o/,-;992 93 
+A2A3 



>4 



EE2^ 

A2A3 92 93 



^„ /,/,+;99293\* /+;99293 

Kfl (V_A2A3 ) V^+A2A3 



A2A3 g2 93 

EE^M^^^f)*^;3:::5f. 

A2A3 92 93 



where we used the notation 

(65) d[23] = [dx]3[d'^kA_]3 3S{x-xi)S^^\k±-ki±). 

Clearly, the TMDs associated to the monopole structures {fl,g\j^,h\) are represented 
by overlaps with no global change of orbital angular momentum /S.(.z = 0, the ones 
associated to the dipole structures {f^r^ ,g1j,,h^'^ ,h^^) involve a change by one unit of 
orbital angular momentum | A£^ | = 1 and the one associated to the quadrupole structure 
{h^^) involves a change by two units of orbital angular momentum \^iz\ = 2. If one 
neglects gauge-field degrees of freedom, the brackets in Eqs. (64d)-(64g) are real, and 
one obtains vanishing T-odd TMDs. 



(^) In the 3Q approach, the spectator system consists of two quarks. It is straightforward to 
generalize the expression for helicity amplitudes to any kind of spectator system, as the latter 
is integrated out. 
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5'2. Overlap representation of the TMDs in the canonical- spin basis. ~ Most of the 
quark models being originally formulated in the instant form, it is more natural to work 
in the canonical-spin basis instead of the light-cone helicity basis. Since we considered a 
frame where the target has no transverse momentum, there is no difference between target 
light-cone and canonical polarizations. Assuming that the quark light-cone helicity and 
canonical spin are connected by the rotation in Eq. (40), the components of the LCWF 
in the canonical-spin basis V'ctict^ct (with ai ='1^1) a-nd in the light-cone helicity basis 
^x X \ (with Ai = ±) are related as follows (^) 



(66) 



V^AiAaA 



2'^3 Z. ^ 

(TiCr2<73 



i'^ , , D^'{'^* D^'{'> D^'{'> 

T'(Ticr2(^3 0"lAi cr2A2 (73X3 



The correspondence between the components in the two polarization bases is given in a 
more explicit form in app. B. Since £)(i/2)t£)(i/2) ~ j fQj- ^]^q spectator quarks, we find 
the explicit overlap representations in canonical-spin basis 



(67a) 


f! 


0-2CT3 <32 93 


; 


(67b) 


qIl 


<72<73 92 93 


5 


(67c) 


n? 


CT2C3 (72 93 


(67d) 


M JIT 


= y'd[23]>;>;23m 

0'2 0'3 92 93 


hu f,/,t;W2 93y a;992 93 


(67e) 


M Bit 


0"2 0'3 92 93 


^R (^n<T2'T3 ) n'T2'T3 _ 


(67f) 


ki_ ul- 
1/ "1 


= y'd[23]>;>;23m 

0"2<T3 92 93 


f^R (^n'T2<T3 ) no-2^3 


(67g) 




^fd[23])^}'me 

0"2 0'3 92 93 


ku /',/;1';992 93y /T;992 93 


(67h) 


2M'^ "it 


(T2CT3 q2q3 


T;99293A ;i;99293 

^0-2 (T3 ) ^-1(72(73 ' 





with t)^^ = i)l + 2m3- biT ^^'^ ^iT = ^1 ^ Im'' ^it- The functions fl'^^jfliT' '^iL^ ^it ^i'^ 
again related to the TMDs g'^j^,g1^rp,h^^,h^^ according to Eqs. (44) and (45). 

5'3. Spherical symmetry. - Wc now discuss how spherical symmetry restricts the form 
of the wave function in the canonical-spin basis. Spherical symmetry requires the wave 



(^) Note that il'aicr2cr3 cannot be identified in general with the usual rest-frame wave function 
^CTiCT2CT3- They have the same spin structure, but not the same momentum dependence. 
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function to be invariant under any rotation, i.e. 



(68) Y. ["(^' ^)U.[ ["(^' '^)]..4 ["(^' ^)].3.^ ["(^' ^)U' <4<.^ = ^^ 



A 



with the SU{2) rotation matrix u{9, (p) given by Eq. (21). In particular, invariance under 
a (tt, 0)-rotation leads to 



A _ /I NA+O-1+0-2+0-3 „/,A 



(69) VI^^_,^_,3 = (-l)^+-^+-=+-3^, 



while invariance under (0, (/))-rotations implies that all components with i!^ 7^ have to 
vanish 



(70) V'Itt = V^ut = V^It; = ^lu "" '^lu = 0- 



Taking into account the constraints (69) and (70) in an arbitrary (0, 0)-rotation, one 
finally gets (J) 



(71) vjt; + ^nr + ^In = 0. 



Again, spherical symmetry implies that the TMDs are either identically zero or pro- 
portional to the unpolarized and polarized amplitudes A'' and B'^, 



(72a) fl = A\ 

(72b) gl^^ cose B\ 

(72c) hl^'-2^B<^, 



(72d) ^/i7=0, 

(72e) ^5?T-si, 

(72f) ^^]''=0, 

cos 6* — 1 



M •'IT 
M 
M 
M 

(72h) ^/^^T^-^^^^^^i^^ 



(72g) |i-/i^^' = -sin0i?^ 



('') Note that spherical symmetry neither restricts the number of non-zero components of the 
wave function in the light-cone helicity basis nor relates them in a simple way, see Table II in 
app. B. 
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(73) 
(74) 



92 93 
92 93 

'd[23]^[(V' 



99293 |2 _|_ |j,T;99293|2 , |,/,T;99293 |2 

tit I "^ l^itt I 



\4 



^//^'l^ + l^llf^^p-l^ji^^^'^p 



\4 



,t^293) ^1^293^(^1^293) ^. 



i; 992 93 



The TMD relations (36)-(38) then follow trivially. 

5'4. SU(6) spin-flavour symmetry. - Many quark models, in addition of being spher- 
ically symmetric, assume also the SU{6) spin-flavour symmetry. As a result, the wave 
function in the canonical spin basis is given by the product of a symmetric momentum 
wave function (f) and a spin-isospin component ^^fj^ira'' 

(75) i^^'lfS^' = Hk , h , h)'^^-f.\'St ■ 

The wave function in the spin-flavour space can be written as 



(76) 
with 
(77) 
(78) 



^zw'^-^W'^c+x'e] 



V2 



[m-nt], r 



V6 



[2;n-nt-m], e^ 



1 

1 
7^ 



[uud — udu] , 
[2duu — udu — uud] 



As a result, writing explicitly the spin and isospin components of the wave function (75), 
we find the following results in the canonical spin basis 



J;t;9l92 93 



(79) 



tn 
m 
in 



uua u 



du dun 



20 -(j) 
-<t> 20 



-(j) 20 



with = 0({fci}) normalized as /[dxja [d^/cj^]3|0p = 1/6. This implies that the unpolar- 
ized and polarized amplitudes A'^ and B'' are simply proportional 



(80) 



A" = 2A'* = I B" = -&B'^ = 12 I d[23] 101^ 



and so the flavour-dependent relation (39) follows trivially with V^ = B'^ /A'^ 
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6. Quark models 

In this section we review different quark models wliicli liave been used for the calcu- 
lation of TMDs. In particular, we summarize the main ingredients of the models and 
discuss whether they satisfy the conditions of sect. 3. In order to facilitate the discussion, 
we sort the quark models in classes defined as follows: 



• 



• 



• 



The light-cone constituent quark model (LCCQM) of Ref. [75] and the light-cone 
quark-diquark model (LCQDM) of Refs. [38, 39, 83] constitute the class of light- 
cone models; 

The covariant parton model of Ref. [37] constitutes its own class; 

The bag model of Refs. [48, 49] and the light-cone version of the chiral quark-soliton 
model (LCxQSM) of Refs. [23, 84, 85] constitute the class of mean-field models; 

The quark-diquark models of Refs. [28, 47, 29, 31, 35] constitute the class of spec- 
tator models. 



We will not discuss the quark-target model of Ref. [51] as it deals with gluons and 
therefore does already not satisfy the first condition of sect. 3. 

6'1. Light- cone models. - The class of light-cone models is characterized by the fact 
that the target state is expanded in the basis of free parton (Fock) states. One usually 
truncates the expansion and considers only the state with the lowest number of partons. 
In the LCCQM, this lowest state consists of three valence quarks, while in the LCQDM 
it consists of a valence quark and a spectator diquark. 

It is well known that light-cone helicity and canonical spin of free partons are simply 
related by the so-called Melosh rotation [86]. Its j = 1/2 and j = 1 representations [87] 
are given by (see app. A for the definition of the spinors and polarization four- vectors) 

r)(i/2)*/M ^ ULc{k,X)u{k,a) 
•^^ ^ ^ 2m 

^J_fV2k++m -kn 

^ ' VN\ kL V2k++m 

D'^2*{k) = -elc{kA)-e{k,a) 

{y/2k++mf -V2{V2k+ + 7n)kji kj 

(82) =—\V2{V2k+ + 77i)kL (V2fc+ + m)'-fci -V2 {V2k+ + m) kn 






k} V2(V2k+ + m)kL (\/2fc+ + 






where m is the parton mass and A^ — (v^ fc+ + rn)^ + k^. The LCWF in the canonical- 
spin basis being identified in these models with the instant-form wave function, it follows 
that \^k'^ — xMq with JvIq = ^^ uji the mass of the Fock state and uji the free energy 
of parton i. Comparing now Eqs. (81) and (82) with Eqs. (40) and (C.9), we obtain 

ion\ 9 m + xMo , . e k± 

( 83 1 cos § = ^= — and sm § = ^-^ . 
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Finally, both the LCCQM and LCQDM consider wave functions with spherical symmetry 
and SU{6) spin-flavour symmetry. In other words, all the conditions of sect. 3 arc satisfied 
in these models, and so are the TMD relations (36)-(39) (^) . 

6'2. Covariant parton model. - The standard quark-parton model (QPM) refers to the 
infinite momentum frame (IMF), where the parton mass can be neglected. The covariant 
parton model is an alternative to the QPM that is not confined to a preferred reference 
frame. Following the standard assumptions of the QPM, the covariant parton model 
describes the target system as a gas of quasi- free partons, i.e. the partons bound inside 
the target behave at the interaction with the external probe (at sufficiently high Q^) as 
free particles having four-momenta on the mass shell. However, since the covariant parton 
model does not refer specifically to the IMF, the parton mass (^) m is not neglected. 
One also assumes explicitly that the parton distributions are spherically symmetric. 

The covariant parton model does not refer explicitly to quark canonical spin or light- 
cone helicity. Instead, it deals with the covariant quark polarization vector. Identifying 
in the Bjorken limit the Lorentz structures of the hadronic tensor with those of the TMD 
correlator, the authors of Ref. [37] found that the TMDs are given in the covariant parton 
model by 



(84a) 

(84b) 
(84c) 
(84d) 
(84e) 
(84f) 

(84g) 



f! 



91l 



(m + xMf + k 
(m + xMf - k 



{dP}, 
{dfci}. 



k± f 
M J - 



-L<? _ 



M J IT 



0, 



M 



glj, = {m + xM)kx_ {dk^}, 



^h^'^ - n 






hll = - (m + xM) fc_L /{dfci}, 



where M is the target mass, {dfc^} and {dfc^} are the measures associated to the distri- 
butions of unpolarized and polarized quarks, respectively. Comparing with Eq. (72), we 



(*) In the derivations of sects. 4 and 5, we tacitly assumed that the rotation connecting light- 
cone and canonical polarizations depends only on the momentum of the parton under consid- 
eration. The Melosh rotation involves the free invariant mass A^o and therefore the momenta 
of all partons in the state. We are in fact not allowed to pull the factor due to the rotation of 
the active quark out of the integral like in Eq. (72). Nevertheless, this technical detail does not 
affect the conclusion about the validity of the relations (36)- (39). Our tacit assumption, which 
does not apply to this specific case, has been introduced to keep the presentation as simple as 
possible. 

(^) Note that the parton mass appearing in the model has to be regarded as an effective mass, 
in the sense that it corresponds to the mass of the free parton behaving at the interaction like 
the actual bound parton. 
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find that 

a 771 + XM n k± 

(85) cos 2 = — 1^^^^=^^^^^ and sm ^ = 



(m + xM) + fc^ Y (™ + 2;^) + fel 

which is nothing else than the Melosh rotation. This is consistent with the fact that the 
active quark is quasi-free in this model. The difference with light-cone models is that the 
physical mass of the target AI is used in the Melosh rotation instead of the free invariant 
mass Mo (^*'). The conditions 1-3 of sect. 3 are therefore satisfied in the covariant 
parton model, and so are the TMD relations (36)-(38). Since this model does not use the 
language of wave functions, the implementation of SU{6) spin-flavour symmetry is more 
delicate and one has to assume that the unpolarized and polarized distributions become 
simply proportional in order to recover the TMD relation (39). 

6'3. Mean-field models. - In mean-field models, the target is considered as made of 
quarks bound by a classical mean field representing the non-perturbative (long-range) 
contribution of the gluon field. Accordingly, the positive-frequency part of the quark field 
appearing in the definition of the correlator (1) is expanded in the basis of the bound-state 
solutions e~'^"* ¥5„(fc,(7) instead of the free Dirac light-cone spinors e~'^'^'^ ULc{k, X)- 
Moreover, one truncates the expansion to the lowest mode ip = ipi with energy E^iev = Ei. 

In these models, the bound-state solution (^(fc, a) is called the quark wave function. 
This object is clearly different from the LCWF introduced in sect. 5. In particular, 
the former is a spinor while the latter is an ordinary scalar function. It is however 
possible to relate them. Since we consider twist-2 Dirac operators F, only the good 
components of the spinors are involved in the quark bilinear Tp{k,a')Tip{k,<7). Using 
mg(A) = P+ULc(fc,A)/\/2i/2fc+ (see app. A), we find that 

Tp{k,a')Tip{k,a) = ip\k,a')P+j"TP+ip{k,a) 
(86) =5]F,V(fc')^A.(fc)^^^^'^'^^''^^^^'^^ 



AA' 



V2k+ 



where we have defined Fxaik) = UQi^)'pik,a-). In agreement with [88, 89] where one 
boosts explicitly the system in the mean-field approximation to the IMF, F\^{k) can be 
interpreted as the quark LCWF with k^ = xM — E\cv The mass M is identified with 
the nuclcon mass Mn in the bag model and with the soliton mass Mn in the LCxQSM. 
The mean field is assumed spherically symmetric in the target rest frame. It follows 
that the lowest quark-state solution in momentum space takes the form 

(87) ^(fc,-)=(^^5(|fc|)jx^, 

with Xa the Pauli spinor. The functions / and g in Eq. (87) represent the s (€ = 0) and p 
{(, = 1) waves of the bound-state solution. On the other hand, the general 3Q LCWF for 
a spin-1/2 target involves usually s-, p- and d-waves. There is no contradiction between 



(^'') In the covariant parton model, only the active parton is considered on-shell. In light-cone 
models, all the partons are on-shell so that the Fock state itself is off-shell Mo 7^ M. 
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these two statements since / and g describe a single quark in the target and therefore do 
not represent partial waves of total angular momentum. Note also that in the language 
of 3Q LCWF, the s-. p- and d-waves refer to components with iz = 0, ±1 and ±2, 
respectively. This is an abuse of language as partial waves should refer to I and not Iz ■ 
The quark LCWF corresponding to (87) is then given by Q-^) 

^ ^ '^^^^'V2\ feff(lfcl) /(|fc|) + fe.g(|fc|)J^/ 

It describes in particular how canonical spin a and light-cone helicity A arc related. 
Comparing with Eq. (40), we find that 



N " VN 

with N = p{\k\) + 2 ^ /(|fc|).g(|fc|) +g'i\k\). The 3Q LCWF written as JlLi Px^-.(h) 
times the standard SU{6) spin-flavour wave function with target polarization A = ^^ ai, 
is then consistent with spherical symmetry in the canonical-spin basis. All the condi- 
tions of sect. 3 being satisfied in mean-field models, the TMD relations (36)-(39) follow 
automatically. 

6'4. Spectator models. - The basic idea of spectator models is to evaluate the quark- 
quark correlator <f> of Eq. (1) by inserting a complete set of intermediate states and then 
truncating this set at tree level to a single on-shell spectator diquark state, i.e. a state 
with the quantum numbers of two quarks. The diquark can be either an isospin singlet 
with spin (scalar diquark) or an isospin triplet with spin 1 (axial- vector diquark). The 
target is then seen as made of an off-shell quark and an on-shell diquark. Spectator 
models differ by their specific choice of target-quark-diquark vertices, polarization four- 
vectors associated with the axial- vector diquark, and form factors which take into account 
in an effective way the composite nature of the target and the spectator diquark. 

As advocated in Ref. [90], the parton distributions can conveniently be computed 
using the language of LCWFs. The scalar quark-diquark LCWF is defined as 

(90) ^^(fc) (X ULcik, X)ysULciP, A) 



with target momentum P = 



p+ ML Oi 



Wc do not need to specify all the factors 



in the definition as we are only interested in the structure of the wave function in the 
light-cone helicity basis. The scalar vertex is of the Yukawa type J^s = 5s(fc^)l with 
gs{k'^) some form factor. Writing down explicitly the components, one finds 

(91) ^tik)o,'-^h\'''' -'\, 



( ^) Replacing ip(k,a) by the free Dirac spinor u{k,a), one recovers the Melosh rotation given 

by Eq. (81) ul,(X)u{k,a) = \/2i/2fc+ D^Y^>*(fc). 
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Note the striking resemblance with the Melosh rotation matrix of Eq. (81). One can 
similarly define a rest-frame scalar quark-diquark wave function as 



(92) *^ (X M(A:,(7)3^,M(P,est,A) = g,{k^)^2M{E + m) Sa^ 

with target momentum Prcst = iJ^I,0). This wave function is obviously spherically 
symmetric (^^). Furthermore, Eq. (91) suggests that the quark light-cone helicity and 
canonical spin are simply related by a Melosh rotation, as if the quark was free [91]. 
In other words, the conditions 1-3 of sect. 3 are satisfied in scalar diquark models with 
Yukawa-like vertex, and so are the TMD relations (36)-(38). 
The axial-vector quark-diquark LCWF is defined as 

(93) i;tjk)^ULc{k,X)elc^{K,Xn)yyLc{P,A). 

The spectator model of Jakob et al. [28] assumes the following structure for the axial- 
vector vertex y^ = ^"^^ ' 7,5 (7^ + "17) ^'^d the following momentum argument for the 
polarization four- vector K ^ P. The motivation for such a choice is to ensure that, in 
the target rest frame, the diquark spin-1 states are purely spatial. Indeed, the rest-frame 
axial- vector quark-diquark wave function reads in this model 

(94) ^to ^uik,a)e*^{K,aD)y^u{P,.,st,A) = ^-^ ^2M{E + m) [e,, ■ cr)^^. 

It satisfies the constraints (C.13) and is therefore spherically symmetric. Writing down 
explicitly the components of the corresponding LCWF, one finds 

(95) +° V3^ ^ ^ ° ^/3^ 

^++ = -V2 V^+o, VZ_ = V2^1o, ^X- = ^-- - 0, 

the other components being given by ij}-\-\ = (—1)^+-*'+-^" (^aa ) ' 'with A, A = ±i 
and Ad = +1,0,-1. Again, one recognizes the characteristic factors of the Melosh 
rotation [91]. Comparing the structure of the components of the LCWF in Eq. (95) 
with the structure of the components of the LCWF given in Table IV after applying 
the constraints of spherical symmetry in the canonical-spin basis (C.13), one concludes 
that only the quark polarization is rotated. This is in agreement with the fact that 
the momentum argument of the polarization four- vector e^ does not have any transverse 
momentum, and so there is no rotation of the diquark polarization. All the conditions 1-3 
of sect. 3 being satisfied in the axial- vector diquark model of Rcf. [28], the TMD relations 



( ) The rest-frame wave function in Eq. (92) is expressed in terms of canonical spin and therefore 
has the same spin structure as the LCWF expressed in the canonical-spin basis. It follows that 
the constraints due to spherical symmetry discussed in app. C'l apply also here. Furthermore, 
the momentum-dependent part of the wave function in the rest frame does not depend on a 
specific direction. 
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(36)- (38) follow automatically. The flavour-dependent relation (39) can be obtained by 
further imposing SU{6) spin-flavour symmetry to the wave function (^^). 

On the contrary, some versions of the spectator model presented by Bacchetta et 
al. in Ref. [47] do not support any TMD relation. We therefore expect that at least 
one of the conditions 1-3 of sect. 3 is not satisfied. These versions are based on the 
axial- vector vertex 3^^ = ^° v- ' 7''75 and involve the diquark momentum K ^ P — k 
in the polarization four-vector. With these choices, it is found that the condition 3 of 
sect. 3 is not fulfilled since the corresponding rest- frame wave function does not satisfy 
the requirements of spherical symmetry 

(96) ^Ln ^ (e-« • '^).A > 

in accordance with the discussions of Refs. [92]-[94] and the comment in Ref. [47] that in 
this approach the partons do not necessarily occupy the lowest-energy available orbital 
(with quantum numbers J^ ~ \ and L^ = 0.) 

7. 5[/(6)-symmetry breaking in a light-cone constituent quark model 

In this section we apply the formalism introduced in the previous sections for the 
calculation of T-even TMDs from LCWFs to a specific model, namely the light-cone 
constituent quark model. For recent quark-model calculation of the T-odd TMDs we 
refer to [47, 35, 44, 42, 43, 95, 33, 32]. First, we specify the LCWF of (66) which 
assumes S't/(6) symmetry for the spin-flavour component of the wave function. Then, 
we add a component of the LCWF which breaks the SU{Q) symmetry. 

In a first step the percentage of S'[/(6)-breaking terms in the total LCWF is left as 
free parameter, and we discuss in general the effects of these terms on the TMD results. 
This calculation can be reproduced numerically using the Mathematica program which 
can be downloaded from Ref. [96] in both the Windows and Linux versions (^^). In 
particular, this program allows one to reproduce the model calculation for i) the T-even 
TMDs; ii) the spin densities in the transverse- momentum space as function of the quark 
and nucleon polarizations; in) the fcj^ dependence of the x moments of the TMDs, with 
a comparison of the model results with the Gaussian Ansatz. 

In a second step, we study different observables, which are particularly sensitive to 
SU{Q) breaking and can be used to fix the free parameter of the model corresponding to 
the percentage of S'C/(6)-breaking terms in the LCWF of the nucleon. Explicit examples 
of this fitting procedure are given in a second Mathematica file [97]. In particular, in 
this program we present different strategies for fitting the available experimental data 
of three observables: i) the inclusive polarized asymmetry Ai] ii) the ratio F2/F2 of 
the neutron to proton unpolarized structure function of DIS; Hi) the nucleon electroweak 
form factors. In these notes we will discuss, as an example, the fit of the double spin 



(^^) The scalar and axial- vector diquarks represent in principle more than just two quarks. For 
this reason, they have a priori different masses, cutoffs, form factors, . . . When we impose SU{6) 
symmetry, we implicitly consider that the quark-diquark picture originates from a 3Q picture. 
The scalar and axial-vector diquarks then just differ by their spin and flavour structures which 
are uniquely determined by the 517(6) symmetry. 

(^*) In order to run the program, you also need the auxiliary functions available at the same 
web address with the name "Functions-TMD.zip". 
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asymmetry Ai in DIS, leaving for a future work a systematic study to select the best 
fitting procedure. 

7'1. Light-cone constituent quark model. - The S't/(6)-symmetric component of the 
LCWF in the LCCQM has the form in Eq. (66), with the SU{2) matrix relating the 
light-cone helicity and the spin given by the Melosh rotation discussed in sect. 6'1 and 
the wave function in the canonical spin basis given in Eq. (75). For the symmetric 
momentum wave function one assumes the following functional form 



(97) 



0({fcj) = 2(2^)= 



1 UJlLU2i03 

Mo a; 12:2X3 



1/2 



N' 



(Af2+/32)7' 



where N' is a normalization factor, and the scale /?, the parameter 7 for the power-law 
behaviour, and the quark mass m are taken from Ref. [98], i.e. f3 = 0.607 GeV, 7 — 3.4 
and m = 0.267 GeV. According to Schlumpf's analysis [99] these values lead to a very 
good description of many baryonic properties. 

We now introduce an admixture of mixed-symmetry component in the nuclcon wave 
function given by 



(98) 



^A;«i92g3 3^ (,Qg^^A;<jig2g3 + sin S ih^''^'^'"' 



93 



where i^^fa^ST ^^ *^*^ 5C/(6)-symmetric component of Eq. (75) and i}^l%\''^'^^ is a mixed- 
symmetry component. The factor cos^ 6 and sin S give the percentages in the nucleon 
wave function of 5C/(6)-symmetric and S'C/(6)-breaking terms, respectively. For the 
S'C/(6)-breaking component we take an appropriate combination of mixed-symmetry spin- 
isospin wave functions with two momentum wave function of mixed symmetry 



(99) 



i'aXtrak} 



[0^ (x"r - x^eO + ^" (x"C^ + x^C)] 



where the spin and isospin wave functions are defined as in (77), and the momentum- 
dependent part is given by 



(100) ^^({k})=N-^^<j>i{k}), 



b^ = NP 



13' 



:0({fcJ)- 



In Eq. (100) N" and N^ are normalization factors and the Jacobi coordinates are defined 
as 



(101) 



1 

71 



(fc2 - fca), /3 : 



fci. 



As a result, writing explicitly the spin and isospin components of the wave function (100), 
we find the following results in the canonical spin basis 



,7,t:9i92?3 
Vcricr2Cr3 



(102) 



m 
m 



uud 



udu 



-2cf>f' 
cj}^ + ^/3(/)" 



dun 



-2(t)^ 0^TV3^ 

-t- V3<^" (p^ - V3(f>° 
- V^<j)" -2(j)^ 
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with the normalizations J[dx]3 [d2fc_L]3|0"P = /[da;]3 [li^ k±]3\(l)P \'^ = 1/6. 

We note that the mixed-syniinetry components of the nucleon wave function cor- 
respond to a state with 1^ = and satisfy the relation (71) derived from rotational 
invariance. Therefore, the LCWF overlap representation of the TMDs keeps the form 
given in Eqs. (72), with the functions A'' and B'^ calculated from the LCWF in Eq. (98). 
In particular, separating the contribution from the 5'C/(6)-symmctric component, the 
S'C/(6)-breaking term and the interference term we can write 

(103) A" = cos^ ,5 Aly^ + 2 cos (5 sin S A«„j + sin^ 6 A«^, 

(104) B« = cos^ <5 Bl + 2 cos 5 sin 6 Sf„t + sin^ 6 B^^. 



A^ and B^ are the contributions calculated in Eq. (80), while the contribution from 
the /S'[/(6)-breaking term reads 

Ad 

(105) Al. = ^^6[{cj^^f + {rn, 

(106) Sb". = 4(0^)^ Bt = -icl>^f + 3{rr. 

Furthermore, the contribution from the interference between the two components is given 
by 

(107) ^L - -< = ^Sint = iBf^.t = 6%/20/. 

We immediately notice that the admixture of mixed-symmetry components leads to 
the breaking of the flavour-dependent relation (39) between the unpolarized and po- 
larized TMDs, while it docs not spoil the remaining relations among polarized TMDs. 
Furthermore, the up and down quark contribution to the TMDs are no longer propor- 
tional. 

7'2. Results for TMDs. - In this section we report a few results for the TMDs in the 
LCCQM which can be explicitly calculated from the program of Rcf. [96]. In particular 
we discuss the spin densities in the transverse-momentum plane given by 



p{ka:,ky, (A, s±), (A, S±)) = - 



(108) + s\S\h, + sl(2FF - kl5^^)Si^hiT 
where, for a generic TMD j, we introduced the first-x moment defined as 

(109) j{kl) = Jdxj{x,kl). 
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As discussed in sect. 2, the unpolarized TMD /i, the helicity TMD giL, and the transver- 
sity TMD hi in Eq. (108) correspond to nionopolc distributions in the momentum space 
for unpolarized, longitudinally and transversely polarized quarks, respectively. They can 
be obtained from the overlap of LCWFs which are diagonal in the orbital angular mo- 
mentum, but probe different transverse momentum and helicity correlations of the quarks 
inside the nuclcon. 



0.4 

^ 0.2 
> 

o 0.0 



-^ 



-0.2 
-0.4 







-0.4 -0.2 0.0 0.2 0.4 

k, [GeV] 



-0.4 -0.2 0.0 0.2 0.4 

k, [GeV] 



Fig. 2. - Density of quarks in the fc^ plane for net transverse polarization of quarks and proton 
in perpendicular directions. The left and right panel shows the results for up and down quarks, 
respectively. 



All the other TMDs require a transfer of orbital angular momentum between the 
initial and final state. The h^rp TMD describes the distortion due to the transverse 
polarizations in perpendicular directions of the quark and the nucleon [100]. In this case, 
the nucleon helicity flips in the direction opposite to the quark helicity, with a mismatch 
of two units for the orbital angular momentum of the initial and final LCWFs. In Fig. 2 
we show the results for both up and down quarks in the case of a S't/(6)-symmetric 
LCWF. 

The results for definite quark and nucleon polarizations in transverse directions are 
obtained by adding the monopole contributions from /i and hi in Eq. (108). In the case of 
quark and nucleon transversely polarized in perpendicular directions {p{kx, ky, s^, Sy) — 
p{kx,ky, Sy, Sx)), one has from Eq. (108) the sum of the monopole contributions from /i 
and hi, and of the quadrupole contribution from h^j,, ignoring the contributions from 
the T-odd TMDs which are vanishing in absence of gluon degrees of freedom. Since 
the quadrupole distortion induced by h^j, is quite small with respect to /i and hi, the 
density for definite polarization differs slightly from a monopole distribution. 

Among the distributions in Eq. (108), the dipole correlations related to gix and 
hj^j^ have characteristic features of intrinsic transverse momentum, since they are the 
only ones which have no analog in the spin densities related to the GPDs in the impact 
parameter space [101, 102, 36]. In particular, giT and /i^*"^ correspond to quark densities 
with specular configurations for the quark and nucleon spin: giT describes longitudinally 
polarized quarks in a transversely polarized nucleon, while hj^j^ gives the distribution of 
transversely polarized quarks in longitudinally polarized nucleon. Therefore, gix requires 
helicity flip of the nucleon which is not compensated by a change of the quark helicity. 
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Fig. 3. - Quark densities in the k± plane for net longitudinal polarization of quarks in a trans- 
versely polarized proton for up (left panel ) and down (right panel) quark. 



and vice-versa h^j^ involves helicity flip of the quarks but is diagonal in the nucleon 
helicity. As a result, in both cases, the LCWFs of the initial and final states differ by 
one unit of orbital angular momentum and the associated distributions have a dipole 
structure. The results in the S'C/(6)-symmetric version of the LCCQM for the densities 
with longitudinally polarized quarks in a transversely polarized proton are shown in 
Fig. 3. The results for definite quark and nucleon polarizations, are obtained by adding 
the monopole contribution of /i (the contribution from the Sivers function is ignored 
because we do not include gauge-field degrees of freedom). As shown in Fig. 4, the 
sideways shift in the positive (negative) x direction for up (down) quark due to the 
dipole term oc XS'^k^jj-gir is sizable and corresponds to an average deformation (fc") = 

55.8 MeV, and (k^) = -27.9 MeV. The dipole distortion ex As'k'^h^^ in the case 
of transversely polarized quarks in a longitudinally polarized proton is equal but with 
opposite sign, since in our model hj^j^ = —gir- These model results are supported from a 



0.4 

^ 0.2 
> 
O 0.0 

-0.2 

-0.4 




-0.4 -0.2 0.0 0.2 0.4 

A:. [GeV] 




-0.4 -0.2 0.0 0.2 0.4 
k, [GeV] 



Fig. 4. - Quark densities in the k± plane for longitudinally polarized quarks in a transversely 
polarized proton for up (left panel ) and down (right panel) quark. 
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recent lattice calculation [79]-[81] which gives, for the density related to giT, (fc") = 67(5) 
MeV, and (k^) = -30(5) MeV. For the density related to h^^, they also find shifts of 
similar magnitude but opposite sign: (fc^) = —60(5) MeV, and (fc^) = 15(5) MeV. 

In the program of Ref. [96] one can also calculate and simultaneously visualize the 
effects of SU{6) breaking for the spin densities, and for the x and k\ dependence of the 
TMDs. For example, we summarize a few results which can be easily verified following 
the calculation in the program: first of all, one finds that the fej^ dependence of the 
TMDs is not Gaussian. In particular, the fcj^ dependence is the same for up and down 
quark within the S'C/(6)-symmetric version of the model, while the admixture of mixed- 
symmetry components breaks this flavour independence. 

7'3. Observables. - One of the observables which is largely sensitive to effects of 
SU{6) breaking is the inclusive double spin asymmetry Ai in DIS with lepton and target 
nucleon longitudinally polarized. It can be expressed in terms of the unpolarizcd f1{x) 
and polarized gl{x) parton distributions as 

which becomes, in the case of proton and neutron target, 

^ ' ' 4/r(x) + /f(x)' ' Affix) + f^ixy 

where we used isospin symmetry in such a way that the structure functions of the neutron 
are related to those of the proton by interchanging u and d flavours. 

We notice that in a S'f7(6)-symmetric model, one has gi{x) = — 4gf (x), which implies 
Ai = at the scale of the model. At higher scale, ^" 7^ due to evolution, but the 
effects remain small. This is illustrated in Fig. 5 where we show the results for the 
proton and neutron at the initial scale of the model (dashed curves) and after leading 
order (LO) evolution to Q^ = 2.5 GeV^. We also see that experimentally A^ (extracted 
by subtracting deuteron and proton data, or from "^He data, modulo nuclear corrections) 
is found clearly non zero, giving a signature of S'C/(6)-breaking effects. Following the 
calculation in the Mathematica program of Ref. [97], we can use this sensitivity to SU{6) 
breaking to fit the parameter S of the model in Eqs. (103)-(104) to the experimental data 
for Ai. In [97] we give different examples for the strategies which can be used in the 
fitting procedure. Here we discuss only one of them which seems the most suitable for 
the LCCQM. In order to decide the criteria of the fit, we should answer the following 
questions. First, in which x-range and with what accuracy is the model applicable? 
Second, how stable are the results under evolution? A related key question emerging not 
only here but in any nonperturbative calculation concerns the scale at which the model 
results for the parton distributions hold. From the point of view of QCD where both 
quark and gluon degrees of freedom contribute, the role of the low-energy quark models 
is to provide initial conditions for the QCD evolution equations. Therefore, we assume 
the existence of a low scale Qq where glue and sea quark contributions are suppressed, 
and the dynamics inside the nucleon is described in terms of three valence (constituent) 
quarks confined by an effective long-range interaction. In fact, glue and sea quark degrees 
of freedom might be thought of at this low scale to be contained in the structure of the 
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Fig. 5. - The inclusive double spin asymmetry Ai in DIS off proton (left panel) and neutron 
(right panel) as function of x. The theoretical curves are obtained with g1{x) and fi{x) from 
the LCCQM with SU{6) symmetry as follows: both function LO to the (Q^) = 2.5 GeV^ of 
the experiments (solid curves) and both at the low scale of the model (dashed curves). The 
experimental data for Af are from Refs. [107]-[110] and for Al from Refs. [111]-[114] 



constituent quarks, which arc massive objects. The actual value of Qq is fixed evolving 
back unpolarized data, until the valence distribution matches the condition that the 
second moment (a;(Qo))vai, i-e- the momentum fraction carried by the valence quarks, 
is equal to one [106]. Following [122], wc start with the LO-valuc of {x{Q'^))vai = 0.35 
at 10 GeV^ [123] and using LO evolution equations we find the matching condition 
(2^(Qo))vai = 1, at Qo = 420 McV. Although there is no rigorous relation between the 
QCD quarks and the constituent quarks, which requires a more fundamental description 
of the transition from soft to hard regimes, this strategy refiects the present state of the 
art for quark model calculations [103]-[105], and has been validated with a fair comparison 
to experiments [23, 45, 103, 106]. 

As can be seen in Fig. 5, the description of the experimental data is reasonable. 
For X > 0.15 the model describes the Ai data within an accuracy of about 30%. The 
description improves in the valence- a; region oi x > 0.25. Since the model contains no 
antiquark- and gluon-degrees of freedom, it is not surprising to observe that it does 
not work at small x. In conclusion, it can be said that the results are weakly scale 
dependent, and the model well catches the main features of the observable Ai in its range 
of applicability, namely in the valence- a; region. Therefore, we choose to fit the parameter 
describing the admixture of 5'L'^(6)-breaking terms in the nucleon wave function using 
the LO results for Ai of proton and neutron evolved at the average (Q^) = 2.5 GeV^ of 
the available experimental data, and restricting ourselves to the valence region x > 0.3. 
The results of the fit are: S = 7.25 degrees with x^ = 52.41 for a total of iV = 28 
experimental points (Xred ~ X^ /{^ ~ 1) = 1-94). The fitted value of S corresponds to 
a percentage of 98.4% (1.6%) for the S'C/(6)-symmetric (-breaking) component of the 
nucleon wave function. The results of the fit are shown in Fig. 6, with the error band 
corresponding to the la region. We see that with a small percentage of SU{6) breaking 
we are able to give an overall good description of the asymmetry for both proton and 
neutron target in the valence region. 

The results of the fit can be also tested with other observablcs. In particular wc 
consider the double spin asymmetrsy A^j^ which can be accessed in SIDIS of hadrons off 
proton target. For a more extensive discussion about the calculation within the LCCQM 
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Fig. 6. - The inclusive double spin asymmetry Ai in DIS off proton (left panel) and neutron (right 
panel) as function of a;. The theoretical curves are obtained from LO evolution to Q^ = 2.5 GeV^ 
of the PDFs calculated in the LCCQM with admixture of SU{6) components in the nucleon 
wave function. The parameter describing the percentage of the S'f/(6)-symmetric and SU{6)- 
breaking components is obtained from the fit to the experimental data, with the result 5 = 7.25. 
The experimental data for A^" are from Refs. [107]-[110] and for AI from Refs. [111]-[114] 



of the azimuthal asymmetries in SIDIS related to T-even TMDs we refer to [45] , while 
for comprehensive reviews of recent and planned experiments to access the TMDs we 
refer to [27, 124, 125]. Assuming the Gaussian Ansatz for the k']_ dependence of the /i 
and gi TMDs, this asymmetry can be written as 



(112) 



A 



LL 



Fll ^ T.a^lxgl{x)Dl{z) 
Fuu Y.a^lxf1{x)Dl{z) 



which reduces to Ai in Eq. (110) if no hadrons are observed in the final state. In Eq. (112), 
Di{z) is the unpolarized fragmentation function describing the hadronization process of 
the struck quark decaying into the detected hadrons, with z the energy fraction taken out 
by the detected hadron. In the calculation we used for Di{z) the LO parametrization 
[115] at Q^ = 2.5 GeV^. In Fig. 7, we show the results for All hi DIS production 
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Fig. 7. - The semi-inclusive double spin asymmetry All in DIS of pions off a proton target as 
function of x. The dashed curves correspond to the LCCQM results in the 5'f/(6)-symmetric 
version, while the solid curves include the contribution from 5'f/(6)-breaking terms, with a 
percentage determined from the fit of the model to the A\ (see text). For the fragmentation 
function Di we used the parametrization at Q^ = 2.5 GeV^ from Ref. [115]. The data are from 
Refs. [116, 117]. 
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Fig. 8. - The electric form factor of the nucleon. The dashed curve shows the prediction within 
the S'[/(6)-symmetric version of the LCCQM. The sohd curve is the LCCQM results with the 
inclusion of 1.6% of S'[/(6)-breaking component in the nucleon wave function. The references 
to the data can be found in [121]. 



of pions off proton target. In particular, the theoretical curves are obtained with the 
model results of gi{x) and fi{x) evolved at LO to Q^ ~ 2.5 GeV^. The dashed curves 
show the results from the 5f/(6)-symmetric version of the LCCQM, and the solid curves 
correspond to the LCCQM predictions with a percentage of 5t/(6)-breaking terms as 
obtained from the fit of Ai. In both cases, we show the predictions only for x > 0.3, 
corresponding to the range where we performed the fit for Ai. We clearly see that the 
inclusion of S't/(6)-breaking terms leads to a quite good agreement with the data in the 
valence region. 

An other observable which is particularly sensitive to SU{6) breaking is the electric 
neutron form factor G^. It is well known that in the non-relativistic SU{6) constituent 
quark model G^ is identically equal to zero. The inclusion of relativistic effects, like in 
the LCCQM, produces non-zero results. However, the predictions remain small and do 
not reproduce the behaviour of the experimental data, as shown by the dashed curve in 
Fig. 8. Corrections can come from the meson-cloud of the nucleon [118, 119, 120]. The 
meson-cloud contribution was recently calculated within the LCCQM in Ref. [120] and 
it was found quite smooth, significant only at Q^ < 0.5 GeV^, but not able to reproduce 
the trend of the data. Only the inclusion of S't/(6)-breaking terms can give a substantial 
improvement. This is shown by the solid curve in Fig. 8, obtained within the LCCQM 
with the percentage of S'[/(6)-breaking terms which was obtained from the fit to the 
asymmetry Ai. Remarkably, the good agreement with the experimental data provides 
an important consistency check for our estimate of the S't/(6)-breaking contribution to 
the nucleon wave function. 



8. Conclusions 



In this work we presented a study of the transverse-momentum dependent parton 
distributions in the framework of quark models. In particular we discussed several quark 
models, introducing the different formalisms for the practical calculation of the TMDs 
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and identifying the common building blocks besides the specific assumptions for modeling 
the quark dynamics. We sorted these models in different classes: the light-cone models, 
the covariant parton model, the mean-field models and the spectator models. Most 
of these models predict relations among the leading-twist T-even TMDs. In particular, 
there are in total four independent relations among the leading-twist T-even TMDs: three 
of them are flavour independent and connect polarized TMDs, while a fourth flavour- 
dependent relation involves both polarized and unpolarized TMDs. Since in QCD the 
eight TMDs are all independent, it is clear that such relations should be traced back 
to some common simplifying assumptions in the models. First of all, it was noticed 
that they break down in models with gauge-field degrees of freedom. Furthermore, most 
quark models are valid at some very low scale and these relations are expected to break 
under QCD evolution to higher scales. Despite these limitations, such relations are 
intriguing because they can provide guidelines for building parametrizations of TMDs to 
be tested with experimental data and can also give useful insights for the understanding 
of the origin of the different spin-orbit correlations of quarks in the nucleon. We have 
shown that these model relations have essentially a geometrical origin, and can be traced 
back to properties of rotational invariance of the system. In particular, we identified 
the conditions which are sufficient for the existence of the flavour-independent relations. 
They are: 

1. the probed quark behaves as if it does not interact directly with the other partons 
(i.e. one works within the standard impulse approximation) and there are no 
explicit gluons; 

2. the quark light-cone and canonical polarizations are related by a rotation with axis 
orthogonal to both the light-cone and quark transverse- momentum directions; 

3. the target has spherical symmetry in the canonical-spin basis. 

For the flavour-dependent relation, one needs a further condition for the spin-flavour 
dependent part of the nucleon wave function. Specifically, it is required 

4. SU{6) spin-flavour symmetry of the wave function. 

On the basis of the above assumptions, we were able to derive the model relations among 
TMDs within two different approaches. 

The first approach is based on the representation of the quark correlator entering the 
definition of TMDs in terms of the polarization amplitudes of the quarks and nucleon. 
Such amplitudes are usually expressed in the basis of light-cone helicity. However, in or- 
der to discuss in a simple way the rotational properties of the system, we introduced the 
representation in the basis of canonical spin. In this framework, we showed that the con- 
ditions 1-3 are sufficient for the existence of all three flavour-independent relations. We 
also showed that a subset of these three relations can be derived relaxing the assumption 
of spherical symmetry and using the less restrictive condition of axial symmetry under a 
rotation around a speciflc direction. 

The second approach is based on the representation of TMDs in terms of quark 
wave functions. In particular, we expressed the TMDs as overlap of light-cone wave 
functions, and we derived the relation with the corresponding representation in terms of 
overlap of wave functions in the canonical-spin basis. After discussing the consequence 
of spherical symmetry on the spin structure of the wave function, we were able to obtain 
an alternative derivation of the relations among polarized T-even TMDs. Finally, for the 
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remaining relation among polarized and unpolarized T-even TMDs, we used the SU{6) 
symmetry for the spin-isospin dependence of the nucleon wave function. 

On the basis of this study, we have shown how and to which extent the conditions 
1-4 are realized in the different classes of quark models. In particular we verified that all 
the models satisfying the TMD relations also satisfy the above conditions, while models 
where the TMD relations do not hold fail with at least one of the above conditions. 

As specific example we discussed in details a light-cone constituent quark model, 
giving the necessary ingredients to perform the calculation of the T-even TMDs. The 
calculation can be reproduced numerically using two Mathematica programs which are 
available on-line. In particular, the first program allows one to reproduce the model 
calculation for i) the T-even TMDs; ii) the spin densities in the transverse-momentum 
space as function of the quark and nucleon polarizations; in) the k^ dependence of the x 
moments of the TMDs, with a comparison of the model results with the Gaussian Ansatz. 
The calculation can be reproduced using the S't/(6)-symmetric version of the light-cone 
constituent quark model as well as introducing the effects of 5[/(6)-symmetry breaking 
in the light-cone wave function which lead to the violation of the relation between un- 
polarized and polarized TMDs. The second Mathematica program allows one to obtain 
predictions for different observables, which are particularly sensitive to SU{6) breaking 
and can be used to fix the free parameter of the model corresponding to the percentage 
of S't/(6)-breaking terms in the LCWF of the nucleon. In particular, in this program we 
present different strategies for fitting the available experimental data of three observables: 
i) the inclusive polarized asymmetry Ai; ii) the ratio F2/F2 of the neutron to proton 
unpolarized structure function of DIS; Hi) the nucleon electroweak form factors. In these 
notes we discussed, as an example, the fit of the double spin asymmetry Ai in DIS. 
Using the results from this fit, we also discussed predictions for the semi- inclusive double 
spin asymmetry All in DIS of pious off proton target and for the electric neutron form 
factor. The good agreement of these observables with the experimental data provides an 
important consistency check for our estimate of the S't/(6)-breaking contribution to the 
nucleon wave function. 

This exercise provides a good example of the practical value of models in phenomeno- 
logical studies, showing how model parameters related to particular assumptions on the 
quark dynamics can be tuned to describe available experimental data, and then can be 
used to learn new information on the partonic structure of the nucleon. 
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Spinors and polarization four-vectors 

We collect in this Appendix the different types of free spinors and polarization vectors. 
The free canonical Dirac spinor u{k, a) and polarization four-vector e^{k, a) are given by 



(A.l) 

(A.2) 



u{k,a) 



( ^/E + rn t 



m ' m{E + m) 



where Xt ~ ( o)' Xi ~ ( i )' ^"^^ ^^'^ polarization three-vectors are e^^j^ = -^ (^1,— i,0) 
for Sz = ±1, and Eq = (0,0,1) for s^ = 0. The free light-cone Dirac spinor ULci^i^) 
and polarization four- vector ££,p(A;, A) are given by 
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with e±± = -j= (=f1, — i). Both types of spinors and polarization four-vectors coincide in 
the rest frame fcrest = ("7., 0) 



(A.5) 
(A.6) 



^(fcrost,^) =^ MLc(fcrcst,Cr) = \/2m 

£^{krcst,cr) =eLc(^i-ost,cr) = (0,e^). 



The "good" light-cone spinors are the simultaneous eigcnstates of the operator 75 and 
the projector P+ = ^ 7 7^ 



(A.7) P+ug(A) = mg(A), 75Mg(A) = Aug(A), mg(A) = —i ) xa, 



and one can write 



(A., 



P+ = ^7iG(A)4(A). 
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Table II. - Decomposition in the canonical-spin basis iptia2CT3 of the components of the 3Q 
LCWF m the light-cone helicity basis ip^ ^ ^ . The components are grouped according to the 
values of total orbital angular momentum iz- 
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Appendix B. 

Components of the 3Q LCWF in the Hght-cone and canonical polarization 
bases 

Based on Eq. (66), Tabic II shows explicitly how the components of the 3Q LCWF 
in light-cone polarization basis are decomposed in the canonical-spin basis. We used 
for convenience the notations Zi = cos ^ i ^i = ^ii sin -^ and r^ = km sin -^ for the 

components of the rotation matrix D^/^, of Eq. (40). For example, from the first row 
of Table II, we have 



(B.l) 



4't++ = ^1^223 V'j-t-'^ + ^1^2/3 ip\^^ + zihzs V"!^.,. 4- /1Z2Z3 ip^^^ 
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It is interesting to note that any single component of the 3Q LCWF in the canonical-spin 
basis contributes to all components in the light-cone helicity basis, and vice-versa. So 
even if one considers that the wave function has only components with i^ = in the 
canonical-spin basis, the components of the wave function in the light-cone helicity basis 
present all the values £z = ^1, 0, +1, +2, the orbital angular momentum being generated 



by the rotation matrices D 



(1/2)* 



Appendix C. 



Connection to a quark-diquark picture 

We show in this Appendix how the 3Q picture can be connected to a quark-diquark 
picture. In the latter, one considers the whole spectator system as an object with the 
quantum numbers of two quarks, namely a diquark. One may also assume that this 
diquark does not contain any internal orbital angular momentum. From a 3Q picture, 
this amounts to set ^2 = ^3 = k]j/2 and mjj = 2m with k]j and mo the light-cone 
momentum and mass of the diquark, and m the mass of a valence quark. 

C'l. Scalar diquark. - The scalar diquark is obtained by coupling the two spectator 
quarks so to form a system with total angular momentum j = 0. The LCWF of the 
scalar quark-diquark system can be written in terms of the 3Q LCWF as follows 
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The total orbital angular momentum of a given component ip^ is given by the expression 
4 = A- A with A,A = ±i. 

The corresponding LCWF in the canonical-spin basis is defined through 

(C.2) ^t = T.^^D^-l''\ 

a 

and can consistently be written as 

(C.3) ^A(^^^^)^_l^J^A^^(^.^^^^)_^A^^(^^^^^)- 

The explicit decomposition of Eq. (C.2) is displayed in Table III. 



Table III. - Decomposition in the canonieal-spm basis i/)J of the components of the scalar quark- 
diquark LCWF in the light-cone helicity basis tpj^ . The components are grouped according to the 
values of total orbital angular momentum £z . 
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Spherical symmetry in the canonical-spin basis reads 

(C.4) EK^»U'K^,'^)]lA'^a'' = 

A'ct' 

and in particular implies 

in agreement with Eqs. (69), (70) and (C.3). 
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(C.5a) 
(C.5b) 



C'2. Axial-vector diquark. - The axial- vector diquark is obtained by coupling the two 
spectator quarks so to form a system with total angular momentum j ~ 1- The LCWF 
of the axial-vector quark-diquark system can be written in terms of the 3Q LCWF as 
follows 
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The total orbital angular momentum of a given component ip^^ is given by the expres- 
sion £z — A — X — Xd with A, A = ±| and Xd = +1,0, —1. 

The corresponding LCWF in the canonical-spin basis is defined through 
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with the rotation for the axial-vector diquark given by 
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Provided that OoikD) = ^(^d/S), we can consistently write the axial-vector quark- 
diquark LCWF in the canonical-spin basis as 
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The explicit decomposition of Eq. (C.7) is displayed in Table IV. 

Table IV. - Decomposition in the canomcal-spin basis 'i'lrcro "f ^'^^ components of the axial- 
vector quark-diquark LCWF in the light-cone helicity basis '4'x\ ■ The components are grouped 
according to the values of total orbital angular mom,entum iz . 
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Spherical symmetry in the canonical-spin basis reads 



where 
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in agreement with Eqs. (69)-(71) and (C.IO). 
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